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INTRODUCTION 

Accurate analysis of stress-strain behavior is of critical importance in the 
evaluation of life capabilities of hot section turbine engine components such as 
turbine blades and vanes. The constitutive equations used in the finite element 
analysis of such components must be capable of modeling a variety of complex 
behavior exhibited at high temperatures by cast superalloys. The classical 
separation of plasticity and creep employed in most of the finite element codes in 
use today is known to be deficient in modeling elevated temperature time dependent 
phenomena. Rate dependent, unified constitutive theories can overcome many of 
these difficulties and may be more suitable for the analysis of the complex 
behavior of high temperature superalloys. However, many aspects of the unified 
theories have not been fully evaluated. There is an urgent need for a 
comprehensive evaluation and further refinement of the capabilities of unified 
constitutive models for analysis of high temperature superalloy behavior. 

OBJECTIVE 

It is the purpose of this contract (NAS3-23927) to thoroughly evaluate the 
unified constitutive theories for application to typical isotropic cast nickel base 
superalloys used for air-cooled turbine blades and vanes. The specific modeling 
aspects evaluated are: uniaxial, monotonic, cyclic, creep, relaxation, multiaxial, 
notch and thermomechanical behavior. Further development of the constitutive 
theories to model thermal history effects, refinement of the material test 
procedures, evaluation of coating effects and verification of the models in an 
alternate material will be accomplished in a follow-on for this base program. 

APPROACH 

The scope of the overall program covers several aspects of the development of 
constitutive models for material behavior. The objectives of the base program is 
being accomplished through a two year combined analytical and experimental 
program. This is divided into several tasks, each task focusing on a specific 
objective. First an extensive literature survey was made to identify possible 
constitutive models for detailed evaluation. Based on the detailed evaluation, two 
models have been selected for implementation into a finite element code. A 
comprehensive uniaxial smooth specimen material test program is defined so as to 
investigate the constitutive behavior patterns of Rene' 80, which is the base 
material. These experimental results are being used for both the determination of 
the material parameters and further evaluation of the predictive capabilities of 
the two models. 


preceding PAGE BLANK NQT FIUWD 


293 



The models will be evaluated for multi axial analysis capabilities, based on 
multiaxial test data. Two types of multiaxial tests are being performed - 
tension-tension type at a notch root using an extended ISG technique at Michigan 
State University (MSU) (Ref. 1) and tension-torsion type on hollow tubes at the 
General Electric Turbine Technology Laboratories. The notch root behavior 
prediction capability of the models will be evaluated based on several benchmark 
notch verification experiments. This part is similar to the work conducted by 
General Electric on the NASA sponsored Benchmark contract (Ref. 2). 

The capability of the constitutive models to analyze the behavior of an actual 
engine component will be verified by performing a finite element analysis of a 
turbine blade tip, similar to that described in Reference 3. 

PROGRESS 

A. MODEL EVALUATION AND SELECTION 

A comprehensive survey of various unified constitutive theories 
available in the literature has been completed. From the 13 models 
surveyed, 5 theories were selected for detailed evaluation. They are the 
models of (1) Bodner et. al., (2) Krieg, Swearengen, Rohde (3) Miller 
(4) Robinson and (5) Walker (References 4, 5, 6, 7, 8). Each of these 
models was programmed as subroutines in a computer program, which performs 
a simple numerical integration of the basic equations. All these models 
involve a nunfber of material parameters. For the purpose of evaluation of 
the theories, constants available for different materials and temperature 
in the published literature were used. Each model was subjected to a 
variety of appropriate loading conditions, so as to evaluate their ability 
to model several basic aspects of high temperature superalloy behavior. 
These include: (1) strain rate sensitivity (2) creep (3) stress relaxation 
(4) history dependence (5) cyclic hardening/softening (6) anelasticity. 

In addition, the models were evaluated in terms of their complexities in 
numerical implementation and material parameter evaluation. 

During the course of this detailed evaluation, several generic features 
of the models have become more evident, such as the roles played by the 
backstress and drag stress. The numerical difficulties special to each of 
the models have also become apparent. Based on this evaluation process, 
two models were selected for further detailed investigation. These were 
(1) the Bodner Model and (2) a generic backstress/drag stress model. In 
the generic backstress/drag stress model, the specific functional forms 
are being chosen based on the behavior observed in the Rene' 80 test 
program. 

B. EXPERIMENTAL PROGRAM 

All specimens used in this program are cast as solid specimens and 
machined to the desired configuration. Because of the desire to determine 
thin wall constitutive relationships applicable to airfoils, the specimens 
are tubular, with approximately .030 in. wall thickness. Tables la and lb 
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show all the monotonic tests and creep tests completed to date. Table Ic 
shows all the isothermal cyclic tests that have been performed. The tests 
cover a range of strain rates from .002 in/in/min to .2 in/in/min and the 
temperature range from 538C (lOOOF) to (982C) 1800F. 

The tests are specially designed to meet the needs of constitutive 
model development. At the same time, efforts were taken to maximize the 
types of data obtained. For example, at the end of the strain rate 
controlled monotonic tests, a stress relaxation test is performed. The 
cyclic tests have an automatic data acquisition system, which is capable 
of getting up to 200 data points for each hysteresis loop. Some examples 
of this can be seen in Fig. 4. 

The tension-tension multiaxial tests are being done by Prof. J.F. 

Martin at MSU. These tests utilize an axi symmetric notched round bar with 
three indentations at the notch root. Both the hoop and axial strains 
will be measured using the interferometric strain gage, similar to that 
used in the benchmark test program (Ref. 2). 

Data reduction procedures also reflect the special needs of 
constitutive model development. For each test, the elastic modulus is 
first determined, based on the initial stress-strain readings. Then the 
inelastic strain is calculated. Since time is recorded at each data 
point, the time derivatives of all measured quantities is calculated. 

Thus stress rate and inelastic strain rate is calculated at each point 
using a 7 point sliding polynomial technique. All the results are stored 
in a computer file which can be directly used as input in material 
parameter evaluation. 

C. EVALUATION OF MATERIAL PARAMETERS 

It has been widely recognized that one of the major sources of 
difficulty in the use of unified constitutive theories is the 
determination of the material parameters. No generalized procedures of 
determining these material parameters are currently available. 

Considerable effort has been made to develop such a method in the present 
contract. 

The approach that is adopted is to develop a computer program which 
directly uses the various test results as input and generates the various 
material parameters as output. The computer program developed is kept as 
flexible as possible, so that different functional forms can be used. 

Such an approach also assures consistency in the treatment of the various 
test data. However, it should also be noted that, while conceptually 
simple, such an approach can be very challenging, mainly due to the 
non-linear equations involved. Such a computer program has been developed 
for a generic back stress-drag stress model. 
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The generic back stress-drag stress model is described by the following 
set of equations for the uniaxial case: 

sgn {a-^) 1 

h = ■ *^1 2 

z = g-| - 92|e^|Z - R 2 3 


In the above, 

. T 

e = Inelastic strain rate 
= Backstress 
z = Drag stress 

R-| and R 2 are static thermal recovery functions. 
fl and gi are the hardening functions. 
f 2 and g 2 are the dynamic recovery functions. 

Equations 1, 2 and 3 are a set of coupled non-linear differential 
equations. The specific forms for the various hardening and recovery 
functions are significantly different for the various models that have 
been published. The approach taken in this project is to choose those 
forms that appear most appropriate for Rene' 80 behavior. To determine 
the various material parameters involved, an iterative approach is used. 

In this, a set of starting assumptions are made which are subsequently 
relaxed. Then successive non-linear optimizations are performed in 
equations 1, 2 and 3 using the experimentally measured quantities as the 
basis. 

Rene' 80 test data at 982C (1800F) has been analyzed in detail using a 
computer program incorporating the procedure described above. Some of the 
notable results are as follows: 

(1) fi = constant, f 2 * constant appears to work reasonably well for this 
case. The constants in Ri have been found using slow strain rate 
monotonic and creep tests. However, the overall contribution of the above 
term seems extremely small, as compared to the hardening and dynamic 
recovery terms. 
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(2) Fig. la shows the results of the iteration procedure after 5 iterations, 
using only the high e monotonic test (.2 in/in/min). It appears that 
parameters determined using the computer program can reproduce the stress 
strain behavior reasonably well. Fig. lb shows the same result, but using 
only the small strain rate monotonic test (.002 in/in/min). The constants 
for these two strain rates are significantly different. 

(3) The monotonic based constants are not able to predict the cyclic 
behavior. Fig. 2a and 2b show the cyclic loop predictions using monotonic 
based constants. Fig. 3a indicates that softening is continuing. In the 
model, the drag stress equation constants control cyclic softening. 

(4) Fig. 3 shows the results of using constants based on all monotonic tests. 
It is seen that these parameters overpredict the high e tests and 
underpredict the low e tests. Thus, although the model appears good for a 
specific strain rate, it does not seem capable of representing the entire 
strain rate spectrum used here (0.002 in/in/min to 0.2 in/in/min). 

(5) Fig. 4a shows the comparison of the test data and model prediction for a 
cyclic test at .2 in/in/min. The result shown is for the 96th cycle. 

The initial hardening shown in the plot is to be disregarded, because the 
prediction was made for only 2 cycles and not the entire 96 cycles. 

Fig. 4b shows similar results as above for the 0.002 in/in^in cyclic 
test. Both Fig. 4a and 4b indicate that the procedure works well for each 
strain rate. However, the material parameters are significantly different 
for the two cases. Here again, the difference is believed to be caused by 
the drag stress equation parameters, as in the monotonic case. This 
points to the limitations of the particular model in representing a wide 
range of strain rate behavior. 

Current work is evaluating the Bodner model, and extending the analysis to 
lower temperatures where less strain rate sensitivity is anticipated. 

(0) FINITE ELEMENT CODE IMPLEMENTATION 

The 2-D finite element code containing Bodner 's constitutive model has 
been completed and tested. The 2-D finite element code utilizes two 
dimensional constant strain triangles and an incremental initial strain 
iteration technique. To facilitate the simulation of arbitrary load 
histories, the load history is partitioned into piecewise linear segments 
with steady state thermal conditions during each segment. In order to 
simplify input, reduce convergence problems and minimize cost, a dynamic 
time stepping procedure is incorporated. The 3-D finite element code 
using 20 noded isoparametric bricks is currently being developed. 

In order to verify the 2-D finite element code with Bodner 's model a 
number of uniaxial test cases were run and compared with published 
results. (References 9-11). In addition, a large two dimensional model 
(Fig. 5a) of the benchmark notch specimen (Reference 2) was constructed 
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and run with three different loading histories and compared with published 
experimental results. An example of these comparisons can be seen in 
Figure 5b. The overall performance of the finite element code with 
Bodner's model was quite good. The cost of running the code is comparable 
to one using a conventional uncoupled plasticity and creep constitutive 
model . 
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Figure 3 
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Figure 5b. Bodner Model Predictions for Benchmark Notch 
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